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In the limit of low particle density, electrons confined to a quantum dot form strongly correlated 
states termed Wigner molecules, in which the Coulomb interaction causes the electrons to become 
highly localized in space. By using an effective model of Hubbard-type to describe these states, we 
investigate how an oscillatory electric field can drive the dynamics of a two-electron Wigner molecule 
held in a square quantum dot. We find that, for certain combinations of frequency and strength of the 
applied field, the tunneling between various charge configurations can be strongly quenched, and we 
relate this phenomenon to the presence of anti-crossings in the Floquet quasi-energy spectrum. We 
further obtain simple analytic expressions for the location of these anti-crossings, which allows the 
effective parameters for a given quantum dot to be directly measured in experiment, and suggests the 
exciting possibility of using ac-fields to control the time evolution of entangled states in mesoscopic 
devices. 



I. INTRODUCTION 

The study of quantum coherent effects in mesoscopic 
systems, such as quantum dots (QDs), is a subject of 
great current interest, both from the theoretical point of 
view, and because of a growing number of possible exper- 
imental applications. One of the most notable of these is 
the swiftly developing field of quantum computation, in 
which the coherent manipulation of entangled quantum 
states is an essential component. Recent experimental 
successes in detecting Rabi oscillations in QD systems 
driven by ac-fieldad have spurred interest in the use of 
intense ac-fields to coherentlit manipulate the time de- 
velopment of electronic statesu. An exciting possibility is 
to make use the jphenomenon of coherent destruction of 
tunneling (CDTp, in which the tunneling dynamics of a 
quantum system become suppressed at certain parame- 
ters of the field. Tuning the driving field thus provides 
a simple mechanism to localize or move charge within 
the QD on a rapid time-scale by destroying or restoring 
the tunneling between regions of the device, so allowing 
ac-fields to be used as "electron tweezers" . 

In this paper we study the use of ac-fields for this pur- 
pose, by investigating the time-dependent behavior of 
electrons confined to a two-dimensional QD with a square 
geometry, under the influence of a strong driving field. 
We use an effective model of Hubbard-type to describe 
the system, which gives a considerable computational ad- 
vantage over standard numerical approaches, and also 
allows us to easily include the important effects of the 
electron correlations produced by the Coulomb interac- 
tion. By integrating the Schrodinger equation in time, 
we find, that for certain values of the field, tunneling be- 
tween different charge configurations within the QD can 
be quenched extremely well. In particular we find that 
tunneling processes parallel to the field can be destroyed 
while transverse tunneling is left unchanged, resulting in 
an effective decoupling between the two halves of the QD. 
We explajin these findings by making use of the Floquet 
approachQ, and show that the points at which tunneling 
is quenched correspond to anti-crossings between Floquet 



quasi-energies, the locations of which can be found accu- 
rately using a perturbational method. This significantly 
clarifies how time-dependent electric fields can affect the 
charge distribution inside a strongly correlated QD, and 
how they can be used for quantum control. 

At the low electron densities typically present in QDs, 
correlations produced by the Coulomb interaction can 
significantly infiuence the electronic structure. Such 
strongly correlated problems are notoriously difficult to 
treat, and the addition of a time-dependent field compli- 
cates the problem even further. When the mean inter- 
electron separation exceeds a certain critical value, how- 
ever, a considerable simplification occurs, as the Coulomb 
interaction dominates the kinetic energy and drives a 
transition to a quasi-crystalline arrangement which min- 
imizes the total electrostatic energy. In analogy to 
the phenomenon of Wigner crystallization in bulk two- 
dimensional systems, such a state is termed a Wigner 
molecule. As the electrons in the Wigner state are 
sharply localized in space, the system can be naturally 
and efficiently discretized by placing lattice points just 
at these spatial locations. A many-particle basis can 
then be constructed by taking Slater determinants of 
single-particle states defined on these lattice sites, from 
which an effective Hamiltonian of Hubbard-type can be 
generated to describe the low-energy dynamics of the 
systemH. A major advantage of this technique over stan- 
dard discretizationO schemes, in which a very large num- 
ber of lattice points is taken to approximate the contin- 
uum limit, is that the dimension of the effective Hamilto- 
nian is much smaller (typically by many orders of magni- 
tude), which permits the investigation of systems which 
would otherwise be prohibitively complex. This approach 
has proven to be extremely successful in treating a vaci- 
ety of static problems, including one-dimensionai|DDsEl, 
two-dimensional QDs with polygonal boundariesQIa, and 
electrons confined to quantum ringaJ. We further develop 
this method in this work by including a time-dependent 
electric field, and study the temporal dynamics of the 
system as it is driven out of equilibrium. 
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FIG. 1: (a) Ground-state charge-density for a two-electron 
square QD. GaAs material parameters are used, and the side- 
length of the QD is 800 nm, placing it in the Wigner regime. 
The dark areas indicate peaks in the charge-density, (b) Lat- 
tice points used for the effective lattice-Hamiltonian. 



II. MODEL AND METHODS 

We consider a system of two electrons confined to a 
square QD with a hard-wall confining potential - a sim- 
ple representation of a two-dimensional semiconductor 



QD. Such a system can be produced by gating a two- 
dimensional electron gas confined at a heterojunction in- 
terface, and by placing a_|gate split into four quadrants 
over the heterostructurcEEl, the potentials at the corners 
of the QD can be individually regulated. In Fig.^ we 
show the ground-state charge-densitss obtained from the 
exact diagonalization of a square QDlI, for device param- 
eters placing it deep in the Wigner molecule regime. It 
can be seen that the charge-density is sharply peaked at 
four points, located close to the vertices of the QD. This 
structure arises from the Coulomb interaction between 
the electrons, which tends to force them apart into diag- 
onally opposite corners of the dot. As there are two such 
diagonal states, degenerate in energy, we can understand 
the form of the ground-state by considering it to be es- 
sentially a superposition of these two states (with a small 
admixture of higher energy states). The four points at 
which the peaks occur define the sites on which the ef- 
fective lattice-Hamiltonian operates, as shown in Fig.|l|b. 

We take an effective lattice-Hamiltonian of the form: 



J 



V niUj 



(1) 



where c\^/cj^ are the creation/annihilation operators for 

an electron of spin a on site i, = c\^c^^, and iii is the 
total charge occupation of site i. The quantity i denotes 
the hopping between adjacent sites, and throughout this 
work we set t and fi equal to one, and measure all en- 
ergies in units of t. V represents the Coulomb repulsion 
between electrons occupying neighboring sites, and U is 
the standard Hubbard-JJ term, giving the energy cost for 
double-occupation of a site. Ei{t) denotes the electric po- 
tential at site i, which in general can have a static and a 
time-dependent component. In experiment, static offsets 
can arise either from deviations of the confining potential 
of the QD from the ideal geometry, or by the deliberate 
application of gating voltages to the corners of the QD. 
Applying corner potentials in this way could be used to 
enhance the stability of the Wigner molecule state, espe- 
cially if the confining potential is softer than the hard- 
wall potential considered here. Corner gates may also 
be used to ensure that the multiplet of states included 
in this effective lattice-model is well-separated from the 
other excited states of the QD system, which will there- 
fore not influence the system's dynamics. In this work, 
however, we do not explicitly consider the effects of static 
gates, and we neglect the influence of small, accidental 
offsets encountered in experiment as we expect them to 
havCpCmly minor effects, and indeed may even stabilize 
CDTB. For convenience, we consider applying an ac- 
field aligned with the x-axis of the QD, which can be 



Ef 



-— cos ujt (2) 



parameterized as: 

E 

Ea = Ed = — cos tot, 

where A, B, C, D label the sites as shown in Fig.|l|b. We 
emphasize that although we have the specific system of 
a semiconductor QD in mind, the effective-Hamiltonian 
we are using can describe a wide range of physical sys- 
tems, including 2x2 arrays of connected QDstJ and the 
quajitum cellular automata systems studied by Lent et 
at3, and our results are thus of general applicability. 

We study the dynamics of the system by placing it 
in a certain initial state, and then integrating this state 
in time under the influence of the effective Hpiiiltonian 
(Q), using a fourth-order Runge-Kutta methodllj. During 
the time-evolution, typically of the order of 50 periods of 
the driving fleld, we measure physical quantities such as 
the particle occupation of the sites ni(t), and also en- 
sure that at all times the unitarity of the wavefunction is 
accurately preserved. As we consider a two-electron sys- 
tem, its eigenstates are symmetric (antisymmetric) under 
particle exchange, corresponding to their singlet (triplet) 
symmetry. The Hamiltonian (|l|) contains no spin-flip 
terms since measurements on semiconductor QDs show 
that—the spin-flip relaxation time is typically extremely 
lon^, and so the singlet and triplet sub-spaces are com- 
pletely decoupled. Thus if the initial state possesses a 
deflnite parity this will be retained throughout its time 
evolution, and we only need to include states of the same 
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FIG. 2: Schematic representation of the two-particle basis 
states for the singlet sub-space of the Hamiltonian. The 
ground state of the QD is approximately a superposition of 
states (1) and (2). 



parity in the basis. We choose to use initial states with 
singlet symmetry, which corresponds to the symmetry of 
the system's ground-state. Simple state counting reveals 
that the singlet sub-space has a dimension of ten, and 
can be spanned by the six states shown schematically in 
Fig.§, together with the four states in which each site is 
doubly-occupied. 

Since the Hamiltonian (|^) is periodic in time, we 
can use the Floquet theorem to write the solutions 
of the time-dependent Schrodinger equation as ip{t) — 
exp(— zejt)^j(i), where ej is called the quasi-energy, and 
(f>j{t) is a function with the same period as the driving 
field, called the Floquet state. This type of expression is 
familiar in the context of solid-state physics, where spa- 
tial periodicity permits an analogous rewriting of the spa- 
tial wavefunction in terms of quasi-momenta and Bloch 
states (Bloch's theorem). The Floquet states, and their 
corresponding quasi-energies, can be obtained from the 
eigenvalue equation: 



(3) 



which, as we show in the Appendix, is the key to ob- 
taining an analytic expression for the quasi-energies by 
means of perturbation theory. 

Using the Floquet states as a basis, the time-evolution 
of a general state driven by the periodic field may be 
written as: 



m = 



Cj e 



(4) 



hence evaluating the quasi-energies provides a simple and 
direct way of investigating the long time-scale behavior 
of the system. In particular, when two quasi-energies are 
close to degenerate the time-scale for tunneling between 
the states becomes-extremely long, producing the phe- 
nomenon of CDTuIIj. 

As the ac- field is aligned with the x-axis of the QD, the 
Hamiltonian (|^) is invariant under the composite parity 
operation x ~x;t — > t + T/2. As a result the Flo- 
quet states can also be classified into parity classes, de- 
pending whether they are odd or even under this parity 
operation. Quasi-energies belonging to different parity 
classes may cross, but if they belong to the same class 
the von Neumann- Wigner theorem forbids this, and con- 
sequently at close approaches they form anti-crossings 
instead. Identifying the presence of crossings and anti- 
crossings in the quasi-energy spectrum provides a neces- 
sary (though not sufficient) condition for CDT to occur. 

The locations of these close approaches between quasi- 
energies piay be found by an analytic approach used by 
Holthauala to treat_aon-interacting electrons, and later 
generalized in Ref.liJ to include interactions. In this 
method the quasi-energies are obtained by first solving 
the Floquet equation ^) in the absence of the tunneling 
terms, and then performing perturbation theory with the 
tunneling terms acting as the "perturbation" . A full de- 
scription of the method is given in the Appendix, and it 
has been shown to work extremely well when the tunnel- 
ing terms are small in comparison to the other energy- 
scales of the problem. 

A convenient numerical method to obtain the quasi- 
energies and Floquet states is to diagonalize the unitary 
time-evolution operator for one period of the driving field 
U{t + T,t). It may be easily shown that the eigenvectors 
of this operator are equal to the Floquet states, and its 
eigenvalues are related to the quasi-energies via = 
exp[— jEj-T]. This method is particularly well-suited to 
our approach, as U(T, 0) can be obtained by propagating 
the unit matrix in time over one period of the field, using 
the Runge-Kutta method described earlier. 



III. RESULTS 



A. Non-interacting case 



which is formally analogous to the standard expansion 
in the eigenvectors of a time-independent Hamiltonian. 
Indeed, in the adiabatic limit, T = 27r/w oo, the 
quasi-energies evolve to the eigenenergies, and the Flo- 
quet states to the eigenstates. An important property 
of expanding in Floquet states (|^) is that it provides a 
extremely valuable separation of time-scales. Although 
the Floquet states explicitly depend on time, they are 
periodic with the same period as the driving field and so 
just influence the dynamics on short time-scales. Con- 
sequently, the long time-scale dynamics of the system 
is essentially determined by just the quasi-energies, and 



We begin our investigation by first considering the sim- 
plest case, that of non-interacting electrons {V ^ 0, U = 
0) . In Fig.^ we show the time-evolution of the system for 
two strengths of electric field at a frequency of w = 8, in 
each case using state (6) (see Fig.|^) as the initial state. 
A consequence of using this initial state is that as the 
ac-field is aligned with the x-axis of the QD, it does not 
break the reflection symmetry between the upper and 
lower halves of the QD, and thus throughout the time 
evolution ua — njj, and ns ~ nc- For the first case, 
when the electric potential E = 30.0, it can be seen that 
the occupation of sites cycles between zero and one, as 
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FIG. 3: Time development of the non-interacting system for 
CO = 8: (a) electric potential, E = 30.0 (b) E = 19.24. Solid 
line indicates the occupation of sites A and D, the dotted line 
the occupation of sites B and C. 



the two electrons perform spatial Rabi oscillations be- 
tween the left side of the system and the right side. The 
picture in the second case, however, for E = 19.24, is 
radically different, with the occupation of sites A and D 
varying little from its initial value of one, and sites B 
and C remaining almost empty throughout the time evo- 
lution. It thus appears that at this second value of E 
the tunneling between the left and right sides has been 
considerably suppressed. 

To confirm that CDT is occurring, we present in Fig.^ 
a comparison of the amplitude of the oscillations of the 
occupation number of site A with the quasi-energy spec- 
trum of the system as a function of the electric potential 
E. It can be clearly seen that for most field strengths 
the charge oscillation has an amplitude of approximately 
one, except at sharply defined minima where it is heav- 
ily quenched. The positions of these minima correspond 
precisely to the locations of exact crossings of the quasi- 
energies. Using the perturbative method described in 
the Appendix reveals that the quasi-energies fall into 
three bands. The central band has quasi-energies of 
e± = zt2Jo{E/u!) (where Jq is the Bessel function of 
the first kind) and cq = 0. The upper band has the 
same quasi-energies, but increased by a constant amount 
of +2, and similarly the quasi-energies of the the lower 
band are decreased by 2. Plotting these quantities in 
Fig.^ demonstrates that the agreement of the perturba- 
tive result with the exact results obtained from the diago- 
nalization of U(T, 0) is extremely good, and corroborates 
our observation of CDT at E — 19.24, as at this point 
E/u! — 2.40 which is the first zero of Jq. This depen- 
dence of CDT on Jo is Ikmiliar from the behavior of the 
driven two-level systemQ'Ej, in which the phenomenon of 
CDT was first noted. 

It is interesting to observe that the quasi-energy spec- 
trum resembles that of non-itttemcting electrons driven 
by an ac-field in a superlatticetj'E^, with the quasi-energy 
crossings corresponding to "miniband collapse" . Indeed 
the lattice model we study can be considered to be a four- 
site chain with periodic boundary conditions in space. 
The difference between the cases arises, however, be- 
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FIG. 4: (a) Quasi-energies of the non-interacting system for 
uj — 8: circles = exact results, lines=perturbative solution, 
(b) Amplitude of oscillation of the occupation of site A. 



cause in the case of a superlattice all intersite tunneling 
processes are suppressed, as the electric field is always 
parallel to the axis of the superlattice. Although our 
system is topologically equivalent to a chain, the fact 
that it has a two-dimensional geometry means that only 
tunneling processes parallel to the field are suppressed 
(A B,C ^ D), while tunneling in directions perpen- 
dicular to the field {A ^ D,B ^ C) is unaffected. 



B. Interacting electrons, no double-occupancy 

We now consider the effect of introducing interactions 
between the electrons, and begin by taking the Hubbard 
U-teriR to be infinitely large - that is, we work in the sub- 
space of states with no double occupation. Our Hilbert 
space is thus six-dimensional, and we use as the basis the 
states shown in Figj^. 

We show the quasi-energy spectrum of this system, 
again for a frequency of w = 8 in Fig.|^. In contrast to 
the non-interacting case, we see that the system presents 
two different regimes of behavior. The first of these is the 
weak-field regime, E < V, a,t which the driving field does 
not dominate the dynamics. In this regime the quasi- 
energy spectrum, and correspondingly the amplitude of 
oscillations, shows little structure, and it is difficult to 
obtain analytical results as the perturbational approach 
is not valid when the tunneling terms are comparable in 
magnitude to the electric field.. 

The second regime occurs at strong field strengths. 
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FIG. 6: Time development of the system for U infinite, V = 
80 and a; = 8: (a) electric potential, E = 100.0 (b) E = 
115.7. Solid line indicates the occupation of sites A and D, 
the dotted line the occupation of sites B and C. 



FIG. 5: (a) Quasi-energies of the system for U infinite, 1^ = 80 
and u = 8: circles = exact results, lines=perturbative solu- 
tion [±2J\q{E / uj)\. (b) Amplitude of oscillation of the occu- 
pation of site A. (c) Detail of quasi-energy spectrum, showing 
an anti-crossing, (d) Detail of amplitude of oscillations. 



E > V, for which the quasi-energy spectrum clearly 
shows a sequence of close approaches. In Fig.||c we 
show an enlargement of one of these close approaches 
which reveals it to be an anti- crossing. Employing the 
perturbation theory demonstrates that the two quasi- 
energies involved in these anti-crossings are described by 
±2Jn{E/tj), where n is equal to V/uj. We may thus 
think of n as signifying the number of photons the sys- 
tem needs to absorb to overcome the Coulomb repulsion 
between the electrons occupying neighboring sites. To ex- 
amine whether these anti-crossings correspond to CDT, 
we follow the same procedure as before, and study the 
amplitude of oscillations in ua as a function of E, when 
the system is initialized in state (6). The results in Fig.||b 
and Fig.^ strongly confirm that tunneling is highly sup- 
pressed at the anti-crossings, and hence that CDT indeed 
occurs. 

We further show in Fig.^ the time-dependent number 
occupation of the four sites at two values of E. In Fig.^ 
E has a value of 100.0, and it can be clearly seen that 
the electrons perform driven spatial Rabi oscillations be- 
tween the left side of the QD and the right. Accordingly 
the occupation of the sites varies continuously between 
zero and one. In Fig.^ we show the result of chang- 
ing the electric potential to a value of E = 115.7, which 
corresponds to the center of the first anti-crossing. In 
dramatic contrast to the previous case, we see that the 
occupation of site A and D only varies slightly from unity, 
while site B and C remain essentially empty throughout 
the time-evolution. Only a small amount of charge can 
transfer per period of the driving field between the left 
and right sides of the system, producing the small spikes 
visible in this figure. The amplitude of these features 
is extremely small, however, indicating that the tunnel- 



ing between left and right sides has been almost totally 
quenched. The efficiency of the quenching depends on 
the value of V which for the case we consider is rela- 
tively high. For smaller values of V, qualitatively the 
same features occur, but the efficiency of the quenching 
is diminished, and the sharpness of the anti-crossings is 
reduced. 



C. Interacting electrons, double-occupancy 
permitted 

We now take the most general case, and consider the 
competition between the U and V terms. Setting U to 
a finite value means that the four doubly-occupied ba- 
sis states are no longer energetically excluded from the 
dynamics, and accordingly we must take the full ten- 
dimensional basis set. 

Although it is difficult to obtain precise estimates for 
the values of parameters of the effective Hamiltonian, it 
is clear that in general U > V. Accordingly we choose 
the parameters U — 160, y = 16 to separate the two 
energy-scales widely for our investigation. As before we 
set the frequency of the ac-field to w = 8, and in Fig.0a 
we show the quasi-energy spectrum obtained by sweep- 
ing over the field strength. It is immediately clear from 
this figure that for electric potentials E < U the form of 
the spectrum is extremely similar to the infinite-f/ case. 
Performing perturbation theory confirms that, as in the 
previous case, the behavior of the quasi-energies is given 
by ±2Jn{E/uj) where n — V/u;. We show in Fig.|^ the 
amplitude of the oscillations of ua when the system is 
initialized in state (6), which demonstrates that at the 
locations of the anti-crossings the tunneling parallel to 
the field is again quenched. 

When the electric potential exceeds the value of U, 
however, new structure appears in the quasi-energy spec- 
trum. A group of four quasi-energies, that for weaker 
fields cluster around zero, become "excited" and make 
a sequence of anti-crossings as the field strength is in- 
creased. The amplitude of these oscillations is compa- 
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FIG. 7: (a) Quasi-energies of the system for U = 160 and 
V — 16, Lu — 8: circles = exact results, lines=perturbative 
solution [±2J2{E/uj)]. (b) Amplitude of oscillation of the 
occupation of site A, with (6) as the initial state. 



FIG. 8: (a) Quasi-energies of the system for U — 160 and 
V — 16, ui — 8: circles — exact results, lines=perturbative 
solution [±2Jig(i7/a')]. (b) Amplitude of oscillation of the 
occupation of site A, with site A doubly-occupied as the initial 
state. 



rable to the amplitude of the two quasi-energies dis- 
cussed above, but it is clear that the two sets of anti- 
crossings are not in phase with each other. Perturba- 
tion theory predicts that these high-field quasi-energies 
are given by zt2Jjn{E/uj), where m — (U — V)/uj, and 
thus these anti-crossings arise when the absorption of m 
photons equates to the electrostatic energy difference be- 
tween the two electrons being on neighboring sites, and 
doubly-occupying one site. This then indicates that this 
structure arises from the coupling of the ac-field to the 
doubly-occupied states. 

To probe this phenomenon, we time-evolve the system 
from an initial state consisting of two electrons occupying 
site A. In Fig^p3 it can be seen that for electric potentials 
weaker than U the amplitude of the oscillations in ua re- 
mains small, and shows little dependence on the field. 
As the potential exceeds U, this picture changes, and the 
ac-field drives large oscillations in ua, and in fact mainly 
forces charge to oscillate between sites A and B. At the 
high-field anti-crossings, however, the tunneling between 
A and B is suppressed, which shuts down this process. 
Instead, the only time-evolution that the system can per- 
form consists of undriven Rabi oscillations between sites 
A and D, perpendicular to the field. As these oscillations 
are undriven they have a much longer time-scale than 
the forced dynamics, and thus during the interval over 
which we evolve the system (approximately 50 periods of 
the driving field) , the occupation of A only changes by a 
small amount, producing the very sharp minima visible 
in Fig.||b, centered on the roots of Jm(E/uj). 

As the tunneling perpendicular to the field is undriven. 



it is straightforward to evaluate the time evolution of 
the initial state, if we assume that the left side of the 
QD is completely decoupled from the right side. The 
occupation of sites A and D is then given by: 

nA{t) = 1 ^ cos Vtjit, riD = 1 — cos 57 /^t (5) 

where Vtn = At^ /{U — V), which for the parameters we 
use, gives a Rabi period of Tr = i-w/Vlii = 226.19. In 
Fig.^ we display the occupations of sites A and D as a 
function of time, for two values of electric potential. At 
the first value, E = 200, tunneling between the left and 
right sides of the QD is not quenched, and accordingly 
the occupation of the two sites varies rapidly between 
zero and two as the electrons are driven by the ac-field 
around the system. The second value, E = 185.8, cor- 
responds to the first high-field anti-crossing. It can be 
clearly seen that the charge oscillates between sites A 
and D, with a frequency close to ^Ir. These Rabi oscilla- 
tions are damped, however, indicating that the isolation 
between the left and right sides of the QD is not perfect. 
In this sense we can regard the two sites B and C as 
providing an environment, causing the quantum system 
composed of sites A and D to slowly decohere in time. 
When the tunneling between the left and right sides of 
the QD is strong, for example at E = 200, this decoher- 
ence occurs very rapidly. By moving to an anti-crossing, 
however, and suppressing the tunneling, the rate of mix- 
ing between the two sides of the QD can be considerably 
reduced, and is just limited by the separation in energy 
between the two quasi-energies. Tuning the parameters 
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FIG. 9: Time development of the U — 160 and V = 16 system 
for u! = 8: (a) electric potential, E = 200.0 (off-resonance) 
(b) E = 185.8 (on-resonance). Thick solid line indicates the 
occupation of site A, the thick dotted line the occupation of 
site D. Dotted lines in (b) show the Rabi oscillations of the 
isolated two-site system, Eq.Bl 



and states in the single-occupation sub-space, make it 
an excellent candidate as a control parameter for manip- 
ulating the dynamics of strongly correlated electrons in 
mesoscopic systems. 

Acknowledgments 

The authors would like to thank John Jefferson for 
helpful discussions and comments. This research was 
supported by the EU through the TMR programme 
"Quantum Transport in the Frequency and Time Do- 
mains" , and by the DGES (Spain) through grant PB96- 
0875. 



APPENDIX A 



of the driving field therefore gives us a simple and control- 
lable way to investigate how a two-electron wavefunction 
can decohere in a QD. 



IV. CONCLUSIONS 

In this paper we have studied the time-dependent be- 
havior of a system of two electrons confined to a square 
QD, driven by an external ac-field. By considering the 
strongly correlated limit of the system, we are able to 
use an effective lattice model of just four sites, which 
arises from the natural discretization presented by the 
system in a Wigner molecule state. We emphasize, how- 
ever, that the form of the Hamiltonian we use is very 
general, and our results are thus applicable to a wide va- 
riety of mesoscopic systems. In our effective model, the 
inter-electron Coulomb interaction is described by two 
parameters, U and V, and the dynamics of the system 
consists essentially of tunneling from corner to corner, 
along the perimeter of the QD. We find that when the 
frequency of the driving field is in resonance with V or 
{U ~V), charge moves freely around the system, except 
at sharply defined field strengths at which tunneling par- 
allel to the field is destroyed - CDT. By using Floquet 
theory we have established that these points correspond 
to the roots of Jm{E/oj), where m is the order of the 
resonance (i.e. muj = or muj = [U — V)). 

When CDT occurs, tunneling parallel to the electric 
field is suppressed, and the system is restricted to per- 
forming undriven Rabi oscillations in the perpendicular 
direction. We can therefore expect a radical modification 
in the frequency and polarization of the EM radiation 
emitted by the QD in this situation. In experiment the 
phenomenon of CDT should thus be readily measurable, 
allowing the measurement of the U and V parameters, 
and hence giving a simple and useful parameterization of 
the QD system. The tunability of this effect, and its abil- 
ity to discriminate between states of double-occupation 



The full spin-dependent wavefunction of a two-electron 
system can be factorized into a spinor and a purely spa- 
tial wavefunction: 

*(?-i,c^i;»-2,cr2) = x(cri,o-2)V'(»^i,'r2), (Al) 

where for states of singlet symmetry the spatial wave- 
function is symmetric under particle exchange (and the 
spinor x is anti-symmetric), while the reverse is true for 
triplet states. In this work we consider just the singlet 
sector of the model, and the two-particle states shown 
schematically in Fig.|^ represent the spatial components 
of the two-electron wavefunctions. These are formed from 
symmetric combinations of single-particle states defined 
on the lattice-points A, S, C, 13, so that, for example: 

|-^^ ^ V'A(?-i)V'c(r2) +V'c(ri)^yt(r2) ^^^^ 
-v/2 

with similar expressions holding for each of the states 
shown. An analogous set of states with triplet symmetry 
can be constructed by taking anti-symmetric combina- 
tions. To complete the singlet basis we must include the 
four states in which a lattice point is doubly occupied, 
holding one spin-up electron and one spin-down, which 
we label as: 

|7)= VA(ri)V'A(r2) |8)= i'B{ri)i^B{r2) 
|9) = il'c{ri)il^c{r2) |10) = ^_D(ri )?/;£) (»^2) 

Our starting point to obtain approximate expressions 
for the quasi-energies is Eq|3[ As the Floquet states are 
periodic functions of time, it ^-useful to work in a Hilbert 
space of T-periodic functionsEj, by defining an appropri- 
ate scalar product function: 

((0aW|0b(t)» = ^ j\Mt')\Mt'))dt' (A3) 

where the single-bracket denotes the usual inner product 
for the spatial component of the wavefunctions. The ad- 
vantage of working in this extended Hilbert space is that 
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the Floquet states are stationary states of the operator 
H{t) — H{t) — i{d/dt), and thus we are able to employ 
standard stationary perturbation theory techniques, and 
avoid the complication of using time-dependent methods. 

We first divide the Hamiltonian (|l|) into two parts: Ht 
which contains all the tunneling terms, and Hj contain- 
ing all the interaction terms (those involving [/, V and 
the electric field). We then find the eigensystem of the 



operator 7ii{t) = Hj — i{d/dt), and employ the tunneling 
Hamiltonian as the perturbation. By making this divi- 
sion, we can expect the perturbative result to be good 
when tunneling is small in comparison to the other en- 
ergy scales of the problem, and conversely, to break down 
in the limit of weak fields. 

In the basis we have chosen Hj is diagonal, with en- 
tries: 



J 



Hi = diag (o, 0, y, y - E{t), V,V + E{t), U + E{t), U - E{t), U - E{t), U + E(t) 



(A4) 



r 



where E{t) — Ecoscot^ and thus finding the eigensystem 
of Tii{t) reduces to the straightforward task of solving 
the ten first-order differential equations: 



=e,0,(t). 



(A5) 



I 



For j = 1 and 2, ( [A5| ) has the trivial solutions (t>jit) = 
e, = 0. The third and fifth components also do not 



have an explicit time dependence, and so have similarly 
simple solutions: (f>j{t) = exp[j(ej — V)t]. Imposing pe- 
riodic boundary conditions sets the value of the quasi- 
energy, requiring {ej ~ V) = mui where m is an integer. 
The remaining solutions are all similar in form to each 
other, and as an example, 0io(i) is given by: 



= exp 



-i{U~eio)t- 



E . 
i — sm ujt 



(A6) 



Imposing periodic boundary conditions on the this solu- 
tion requires {U — e) = noj, where n is an integer. 

The eigenvalues of Tii are thus (with a two-fold de- 
generacy), V mod (jU (with a four- fold degeneracy), and 
U mod Lo (also with a four- fold degeneracy). These 



represent the zeroth-order approximations to the quasi- 
energies in the perturbational expansion. It can readily 
be shown that their associated eigenvectors (f)j [t) form an 
orthonormal basis set, and, by using standard degenerate 
perturbation theory, the quasi-energies can be obtained 
to first-order by finding the eigenvalues of the perturbing 
operator {{(f)i\Ht\(l)j)) . By using the identity: 



exp [—i(3smujt\ = 



E 



'/m(/9) exp [— iTOo^f] , (A7) 



to rewrite the ei genv ectors which have the form (A6), the 



scalar products (A3) can be evaluated straightforwardly, 
allowing the matrix elements of the perturbing operator 
to be obtained, and the operator to be subsequently di- 
agonalized. 

For the general case, with V = mcu and U = nuj, 
the first-order approximation to the quasi-energies can 
be readily shown to be (with a four- fold degeneracy), 
zt2Jn{E/u}), and ■iz2Jn-m{E / uo) (with each state being 
two- f old de gener ate). The specific cases treated in Sec- 
tions III ^ and |III B| can be treated in a similar way, 
giving the perturbative solutions quoted there. 
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